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Abstract 



Recent Super-Kamiokande data on the atmospheric neutrino anomaly are 
used to test various mechanisms for neutrino oscillations. It is found that 
the current atmospheric neutrino data alone cannot rule out any particular 
mechanism. Future long-baseline experiments should play an important role 
in identifying the underlying neutrino oscillation mechanism. 







The atmospheric neutrino anomaly [|IJ has been confirmed by the Super-Kamiokande 
experiment 0. The observed up-down asymmetries of the detected muons indicate that, 
with the three known neutrino flavours, the most likely solution to this anomaly is — > v T 
oscillations, although significant mixing with v e cannot yet be excluded [[|. 

The phenomenon of neutrino flavour oscillations was first proposed as a consequence 
of nondegenerate neutrino masses f4|. Although this mass mixing of weak eigenstates is 
the most likely mechanism for the observed atmospheric neutrino oscillations, other possible 
mechanisms have been proposed and, as stressed in Ref. it is important to let experiments 
rather than our theoretical prejudice determine which is the correct mechanism. It is in this 
spirit that we undertake the present study. 

The alternative neutrino oscillation mechanisms considered below share a common fea- 
ture: each requires the existence of an interaction (other than the neutrino mass terms) 
that can mix neutrino flavours. In this paper we shall focus on such interactions which are 
mediated by a scalar (J = 0), a vector (J = 1), or a tensor (J = 2) field. Assuming a two- 
neutrino mixing scheme, the — > v T transition probability for each of these possibilities 
may be parametrized as || 



where E v is the neutrino's energy, L is the neutrino's path length (i.e., the distance from 
where the neutrino is produced to the detector), 5j is a parameter specific to the neutrino 
oscillation mechanism, and 8j is the corresponding mixing angle. The subscript J in 6j and 
8 j is simply a label for the different mechanisms, and does not imply that the values of these 
parameters depend on the value of J. Note the distinct energy dependence of the oscillation 
probability for each value of J. This is the key for determining which neutrino oscillation 
mechanism is at work. 

Proposed mechanisms for the J = case include the mass mixing mechanism ||] for 
which the parameter 6q in Eq.(Q) represents the mixing angle between the neutrino flavour 
eigenstates and the mass eigenstates, and the parameter 5q is given by 




(1) 
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. , . Am 2 ml — m? 

5 (mass) = — = 2 4 1 (2) 

where rrij are the masses of the neutrino mass eigenstates. Another possibility for the J = 
case is neutrino oscillations induced by nonuniversal couplings of the neutrinos to a massless 
string dilaton for which 5q is given by 

<5 (dilaton) = - Am 2 — 20(a 2 "^2 — aim 2 ) . (3) 

Here at are the coupling strengths of different neutrino gravitational eigenstates to the 
dilaton field and it has been assumed that the neutrino mass eigenstates and the gravitational 
eigenstates are the same. The angle 9q for this case is therefore the mixing angle between the 
flavour eigenstates and the gravitational eigenstates. Note that the nonuniversal neutrino- 
dilaton couplings constitute a violation of the principle of equivalence. The parameter (f) in 
Eq.([3]) denotes the local Newtonian gravitational potential. It is assumed to be constant 
over the neutrino path length since the dominant contribution to <fi appears to come from 
the great attractor || which has been estimated to be ||: \<fi\ ~ 3 x 10 -5 . With the constant 
4> approximation, the energy dependence in the oscillation probability is the same for this 
mechanism as it is for the mass mixing mechanism. Neutrino oscillation experiments alone 
will not be able to distinguish these two possibilities. For the purposes of this paper, we shall 
refer to them collectively as the scalar mechanism. Note, however, that the mass mixing 
mechanism requires the neutrino masses to be nondegenerate whereas the dilaton induced 
oscillations can take place even if the masses are degenerate. We also observe that, if the 
neutrino-dilaton couplings are universal, i.e., if ai = «2, and if <5o(mass) ^ 0, <5o(dilaton) 
and <5o( m ass) differ only by a multiplicative constant. 

As an example of the J = 2 case, consider neutrino oscillations induced by the equivalence 
principle violation that results from nonuniversal couplings of neutrinos to gravity ||3_ 0f| . In 
this case 62 is given by f| 

5 2 (VEP) = |0|A 7 , (4) 
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where A7 measures the degree of violation of the equivalence principle and, as in Eq.(||), 
</> denotes the essentially constant local gravitational potential. With a constant <f), this 



mechanism is phenomenologically identical [O] to the velocity oscillation mechanism that 



arises from a breakdown of Lorentz invariance [12] for which 62 assumes the form 



5 2 (velocity) = (5) 

where Av — 1*2 — i>i is the difference between the velocities of two neutrino velocity eigen- 
states. We shall refer to these two mechanisms collectively as the tensor mechanism. In the 
former case, the angle 82 corresponds to the mixing angle between the flavour eigenstates 
and the gravitational eigenstates, whereas in the latter case 82 corresponds to the mixing 
angle between the flavour eigenstates and the velocity eigenstates. It should be emphasized 
that the tensor mechanism can lead to neutrino oscillations even if neutrinos are massless 
(or degenerate). These mechanisms were proposed not so much as a competing mechanism 
for the mass mixing mechanism, but to point out that neutrino oscillation experiments could 
be used as high-precision tests of the symmetry principles fundamental to the theories of 
general and special relativity. 

We include the vector mechanism (the J = 1 case in Eq.(|l|)) in our phenomenological 
study. Although there are models which can yield energy-independent neutrino oscillations 
in matter ||TB|^To]|, we are not aware of any such model that can explain the atmospheric 
neutrino anomaly. The aim of our study is not to test any specific model, but rather to 
check if an energy-independent oscillation mechanism is compatible with the atmospheric 
neutrino data. 

To test the three classes of neutrino oscillation mechanisms, we compare the measured 
values of 

data / n \ 

= {NjN e )\ MC W 
and the up-down asymmetry parameters 



y ,_ (A^VA^)|data a _ 



with the corresponding predictions of these mechanisms, assuming — > v T transitions. 
Here N e and are the number of e-like and /x-like events, respectively. N~ v and N£ v are 
the number of a-like events produced in the detector with zenith angle cos < —rj and 
cos0 > 77, respectively. Note that is defined to be negative for upward going directions 
and we have chosen rj = 0.2 in our analysis. The calculational method is identical to that 



described in flllj and will not be reproduced here for environmental reasons. 

Our results for the i?'s and K's are displayed in Figures 1 and 2 together with the 
Super-Kamiokande results 0, 

i?(sub - GeV) = 0.63 ± 0.03 ± 0.05, 
i?(multi - GeV) = 0.65 ± 0.05 ± 0.08, 
F M °- 2 (sub - GeV) = 0.79 ± 0.05, 
yj- 2 (multi - GeV) = 0.56 ± 0.06. (8) 

The experimental results we use correspond to 535 live days of running. For completeness 
we mention that the Super-Kamiokande results for the e-like up-down asymmetries are 

V e °' 2 (sub - GeV) - 1.14 ± 0.07, 
V e a2 (multi - GeV) = 0.92 ± 0.12. (9) 

Finally, note that only statistical errors are given for the up-down asymmetries since they 
should be much larger than possible systematic errors at the moment. We have not shown 
any figure for Y e - it is expected to be about 1 because —>■ v T oscillations have almost no 
effect on the expected number of electron events. 

From Figures 1 and 2 we see that all three mechanisms can fit the data for a range of 
the parameter 8j. The combined \ 2 fit to V M ,Y" e and R yield the allowed regions for the 
neutrino oscillation parameters shown in Figure 3. The x 2 function is defined to be 

xl tm = x 2 (R) + x 2 (Y) } (10) 

where 



x 2 (R) = E 
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The sum is over the sub-GeV and multi-GeV data samples. The measured Super- 
Kamiokande values and errors are denoted by the superscript "SK" and the theoretical 
predictions for the corresponding quantities are labelled by the superscript "th" . The mini- 
mal x 2 f° r t ne scalar, vector, and tensor mechanisms is 5.4, 5.7, and 10.4, respectively, for 
4 degrees of freedom. All three mechanisms are consistent with the data, but the tensor 
mechanism has the worst \ 2 ■ 

For check we have performed a second \ 2 analysis using the SK \ 2 function, 

2 iq2 

2 _ a , P , 2 i 2 (1 Q\ 

X 2 2 ' Xsub-GcV ' Xmulti-GeV) \ 16 ) 



where 
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and similarly for Xmuiti-GeV- H ere a an d (3 are the absolute and relative normalizations 
of the neutrino flux, y^ are theoretical predictions for each zenith angle bin, are the 
experimental data, and we assume the uncertainties are given by a a = oo and up = 0.12 
(we assumed the larger value of ap given in the third reference of 0). It is understood that 
y% (/z-like events) are multiplied by (1 + a)(l + (3/2) while y e a (e-like events) are multiplied 
by (1 + a)(l — (3/2), and \ 2 is optimized with respect to a and (3. With this alternative \ 2 i 
we also found that all three cases (J = 0, 1, 2) gave an acceptable fit to the data. 

The basic reason that the Super-Kamiokande atmospheric neutrino data cannot distin- 
guish these rather different oscillation probabilities is geometrical. Neutrinos from above 
typically travel quite short distances 50 km) whereas neutrinos going up through the 
Earth travel quite long distances (^ 5000 km). The atmospheric neutrino data can be 



explained by assuming that neutrinos from above do not have time to oscillate while neutri- 
nos travelling through the Earth experience averaged oscillations. In fact the data suggest 
that this occurs for both the sub-GeV and multi-GeV energy range (which is, very roughly, 
0.3 GeV ^ E v < 6 GeV). Thus, while atmospheric neutrino experiments can sensitively test 
for oscillations through the zenith angle dependence (or up-down asymmetry) they are really 
only sensitive to averaged oscillations and are consequently not sensitive to the oscillation 
mechanism. The explicit energy dependence of the oscillation probability is only important 
in the region where the oscillations are significant and are not averaged. The proposed 
long-baseline experiments should be much more sensitive to the oscillation mechanism. The 
reason is that the neutrino path length is fixed and at a moderate distance (typically about 
250 - 600 km). Provided that not all of the oscillations are averaged, and provided that 
at least some of the neutrinos have time to oscillate, then a comparison of the detected 
energy distribution with the expected energy distribution should distinguish the oscillation 
mechanisms. 

In summary, we have examined three distinct oscillation mechanisms (J = 0,1,2 in 
Eq.([TJ)) and compared these with the current Super-Kamiokande atmospheric neutrino data. 
We have found that the current data are quite insensitive to the oscillation mechanism 
responsible for the atmospheric neutrino anomaly. To ultimately determine the underlying 
neutrino oscillation mechanism, it is important to do a more controlled experiment. This 
should be possible in the near future with the advent of long-baseline experiments. 
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FIGURE CAPTIONS 



Figure 1: R as a function of 5j (see Eq. ([!])) for the tensor (bold solid line), vector (solid 
line), and scalar (dashed line) mechanisms. Maximal mixing has been assumed in each case. 
The dimensionless variable X is equal to 82 GeV-km, 5i km, and -fy GeV _1 -km, respectively 
for the tensor, vector, and scalar mechanisms (the choice of X for the mass mixing mechanism 
corresponds to the familiar Am 2 /eV 2 ). The horizontal dashed lines correspond to Super- 
Kamiokande's measured R value ± la error. Figure la is for the sub-GeV neutrinos and 
Figure lb is for the multi-Gev neutrinos. 

Figure 2: Same as Figure 1 except Y~ M (see text) is plotted instead of R. 

Figure 3: The allowed region for the mixing parameter 8j and the oscillation parameter 5j 
obtained with the \ 2 function defined in Eq.([H]). Figure 3a, 3b, and 3c are for the scalar, 
vector, and tensor mechanisms, respectively. 
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Figure 1a 




Figure 1b 
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Fig.3b 
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